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>»! Abstract 

^ . The Impulse Compton Profiles (CP's) J{q) and the Ip^)— expectation 

\o ^ 

^r , values for some inert gas atoms(He-Kr) are computed and compared within 

^^ , the Harbola-Sahni (HS), Hartree-Fock (HP) theories and a Self Interaction 

G^ , Corrected (SIC) density functional model. The Compton profiles for excited 

O , states of Helium atom are also calculated. While the calculated CP's are 

^^! found to generally agree, they differ slightly from one another for small values 
of the Compton parameter q and are in good agreement for large q values. 

k> i The <p"> expectation values within the three theories are also found to be 

H ; 

5t , comparable. The HS formalism is seen to mimic HP reasonably well in the 
momentum space, establishing the logical consistency of the former. 



In the phenomenon of Compton scattering, the Compton cross section (of high energy 
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X-ray or 7-ray photon inelastically scattered by electrons in matter) has a direct bearing on 
the electron-momentum density. In the so-termed impulse approximation |I| the Compton 
cross section is proportional to an experimentally observable quantity, viz. the Compton 
profile(CP), related to the electron momentum density, vide: 

/oo rco 
/ l{Px,Py,q)dpxdpy. (1) 

The momentum density 7(jj) is the diagonal {i.e. p' = p) part of the full, reduced first 
order momentum space density matrix T^l^{p\p') connected in turn, to the many-electron 
momentum-space wave function [^ ^{pi,p2, . . . ,Pn) by 

mom\P\P ) 

= N J^*{p,P2, . . . ,pn)<^{p',P2, . . .,PN)d^P2 ■ --d^PN, (2) 

where a sum over spin may also be included. 

Within the independent-electron approximation such as the Hartree-Fock (HF) theory, 
^momip\p') takes the form 

r^2m{p\p') = j:mimip'), 

i 

with i sweeping through the "occupied" states with the occupancies fi and (pi{p), the 
momentum-space orbital, being the Fourier-transform of the coordinate-space orbital ipilf) 
related through (Hartree atomic units employed throughout): 

0.(p1=(^/e^^>.(r1rfV. (3) 

The <p" >- moments are defined in terms of electron momentum density (EMD) distribution 

by: 

/■oo 

<p" >= 47r / p"+2 l{p)dp, - 2 < n < 4, (4) 

Jo 

where ■y{p) = ^ / '-f{p)dQp is the spherically averaged EMD, in turn leading to the spherically 
averaged impulse Compton profile J{q) = 27i J^^ ^{p) p dp. These <p^> expectation values 

2 



sample the interior as well as exterior regions of the EMD and are also related to atomic 
properties. The < 1/p > moment is twice the peak value of the impulse profile J(0); the 
<p^> moment is twice the kinetic energy (= —Etotai, by the virial theorem) while the <p> 
moment is empirically found to be almost proportional to the exact Hartree-Fock exchange 
energy [g. 

The recent Harbola-Sahni approach [^,^ to the atomic structure calculations proffers an 
attractive alternative to the conventional Hartree-Fock description. The HS approach has 
been proven to be successful in giving the total energies 0^ and co-ordinate space properties 
PJToD practically of Hartree-Fock quality. In addition, coupled with local correlation, it also 
describes the excited states of atoms quite accurately [|lT|. This success of the HS formalism 
prompts one for its critical appraisal in the momentum space through Compton profiles 
and and the < p" > expectation values. We also compare these with the corresponding 
quantities calculated within the HF, the HS and the self interaction corrected (SIC) local 
density functional theories. This study is aimed at bringing out how these p— space quantities 
calculated using the local and orbital-independent prescription of HS compare with those 
of the HF theory, which employs a nonlocal potential (in its exact exchange description), 
and of SIC theory in which the effective potential, although local, turns out to be orbital- 
dependent. In the following we first briefly describe the HF, HS and SIC theories to highlight 
the differences among them. 

All the three theories HF, HS and SIC are independent particle theories in which the 
electron orbitals are obtained by solving the equations (Hartree atomic units i.e. h = \e\ = 
m = 1 are used throughout herein), viz. 

-— + ^^i?(r) + ^;,(f)JV'^ = e^^^; i = l,2,....N, (5) 



where Vnir) = Vnudear{r) + J fr-^rd^r' is the Hartree potential and v^ the exchange potential 



r—r\ 

Here p{r) is the electronic density given in terms of orbitals ipil-f) as p{r) = J2i /i|V'i(^)P- The 
differences in HF, HS and SIC precisely lies in the manner in which the exchange potential is 
prescribed in them. In HF which is the exact theory at the "exchange-only" level, as noted 



above, the potential Vx is nonlocal with its action on ipi{r) is given by 

Vx[r)i^Ar) = 2^ / ^^ — ^ d-'r . (6) 

spin i— spin J 

On the other hand, the exchange potential in the exchange-only HS theory is local and 
is prescribed as the work done in moving an electron in the field of its Fermi hole f^: 



Vx{^=Wx{^ = - / Sx-dl, (7) 

where 

^-(^"')= I T^^A^ - ^') d'r\ (8) 

J \ J' — ^' 

is the exchange "electric field" due to the Fermi hole (or the "Exchange hole") px{r,r'). 

In the SIC theory the exchange potential is calculated within the local-density approxi- 
mation (LDA) which is then ad-hodj corrected for its self-interaction on an orbital-by-orbital 
basis ||12[ . Thus the orbital dependent SIC exchange potential is given as 

4 vr '-J \ f — f' \ 4 vr -' 

where Pi(r) =| ipii"^ 1^ is the orbital density. 

In the context of HS potential, it may be noted that as recently established by Holas and 



March [|I^, the Harbola-Sahni exchange-correlation potential can also be calculated from 
the exact second-order density matrix by employing the differential virial theorem. The 
results are consistent with, and provide the mathematical proof of the formalism proposed 
by HS. In addition, it also spells out how the kinetic energy term missing p,|^ from the 
HS potential arises from the differences in the exact kinetic energy density tensor and its 



Slater-Kohn-Sham orbitals counterpart [|1^]. Similar analysis carried out within the Hartree- 
Fock theory [|T^ reveal that the difference between the exact KS exchange-only potential 
and the HS potential is traced back to the differences in the kinetic energy density tensors 
of the HF theory and its local counterpart (as such, this difference is indeed only marginal) 
TM. Further, within a local prescription, it is not clear as to how one incorporates the 



kinetic-energy effects directly in a self-consistent-field (SCF) scheme. 



In this work the spherically averaged Compton Profiles J{q) within the HF theory are 
computed using the near Hartree-Fock(NHF) quality wavefunctions that employ the STO 



(Slater- Type-orbital) bases tabulated by Clementi and Roetti |^ and the values of the 



<p"> moments are from Ref. [ITB]. On the other hand, the orbitals with effective potential 



Wx are obtained by a modified Herman-Skillman code |T^. The calculated Compton profiles 



along with their available experimental [^,^ and accurate theoretical counterparts |2^ are 



tabulated in tables I-IV for inert atomic systems He-Kr while the moments are displayed 
in tables V-VII. The Helium atom in its ground state has a single orbital, hence all the 
CP's within these three exchange-only theories practically coincide (cf. Table I). For Ne, 
Ar and Kr, it is evident from the Tables II-V that for low q— values the CP's differ from 
each other appreciably. For higher q— values these theoretical (HF, HS, SIC) profiles are in 
better agreement with one another. In the low q— region (g < O.Sa-u), the SIC profiles are 
seen to be the largest and the HS profiles are the smallest in magnitude among the three 
theories. Beyond q ~ O.Sa.u. the three profiles cross each other and are in good agreement 
in the asymptotic region. It is observed further that the "experimental" J(g) is fairly well 
estimated by the "exchange only" theories. It is to be noted that an accurate theoretical 
J{q) computation beyond HF, viz. the configuration interaction calculation for A^e, due 
to Tripathi et al. |2^ while improving upon the HF-CP still slightly underestimates the 
experimental J(0) but overestimates the intermediate profile. The higher values of J{q) in 
the SIC theory indicate that the momentum density is localized near the origin p = in 
the SIC formalism. This can also be seen from the of <p^^> and <p^^ > values which, as 
pointed out above, sample the interior region of the EMD. The higher the momentum density 
near the origin the greater are the values of these moments. That these moments have largest 
values in the SIC formalism may be qualitatively explained as follows: The region near the 
origin in the momentum space by Fourier- reciprocity, corresponds to the asymptotic region 
in the position space. In the density functional theory (DFT) the asymptotic decay of the co- 



ordinate space electron density goes as ~ exp{—2J2\emax\ r), where emax is the eigenvalue 
of the highest occupied orbital |2l|. The highest occupied orbital energy eigenvalues for 



these systems in the SIC formahsin are smaller in magnitude compared to their HF and 
HS counterparts. Consequently, the coordinate-space electron density decays slowly in SIC 
than in the HS and HF theories, leading therefore to higher values of <p~^> and <p'^ > 
moments. The larger values of the HF profiles than HS profiles near q = can also be 
explained similarly. In the HF theory (unlike in DFT), all the orbitals decay with the 



same exponent (~ exp{—J2\e']^f!'\ r)) asymptotically (e^^^, here is the highest occupied HF 
orbital energy eigenvalue) [^] which by reciprocity, reflects in the slower decay of '-/{p) in 
the small \p\, resulting therefore in slightly larger values of < 1/p^ >, <l/p> moments and 
Jhf{(i) (near g = 0) as compared to their HS counterparts. Amongst the various moments, 
the agreement among these theories is the best for the <p^ > moments. The HF and HS 
values of this moment are very close, agreeing up to four significant figures in case of Ne 
and Ar, and up to three significant figures for Kr. This agreement is however not surprising 
since this moment is essentially the double of the negative of the total energy (by the virial 
theorem) and the HS theory is known to produce the total energies which are practically 
equivalent to their HF counterparts. Further, the <p> and <p^> moments are also found 
to be comparable in the three theories. The < p"^ > moments within HS agree with their 
corresponding HF and/or SIC values. Thus, the HS theory with its local prescription for 
the exchange potential seen to mimic the Hartree-Fock formalism reasonably well in the 
momentum space. Our study on the detailed structure of the radially projected first order 
reduced density matrix 0] also supports this similarity between the HS and HF density 
matrices: striking similarities are observed in the structure (contours) of reduced first order 
density matrix in the momentum space. 

The HS theory also offers a simpler description of the bound excited states in comparison 
with the HF description. This is because the HS formalism is not based on the variational 
principle, but rather on the physical effect of the Pauli and Coulomb correlations; which has 
prompted us to compute the CPs for the excited states of helium. Of course, the numerical 
HF approach is also perfectly suitable for excited states with the orbitals identified from the 
number of radial nodes (= n — / — 1) of a given radial part Rni{r) of the orbital. However, 
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the appeal of the HS approach is that it is simple to implement than the numerical HF 
scheme, yielding results that are practically equivalent to the latter. The CP's calculated 
for various excited states of Helium are presented in Table IX. The excited states of helium 
atom will have diffused electron distribution in the position space and also will have higher 
total energy. Consequently, the excited state CP's will be more compact or localized in the 
small q— region, as is evident from Table IX. 

One naively expects that the HS Compton profiles may be improved by adding an accu- 
rate local correlation to its effective potential. It is observed, however, that addition of an 
ad-hoc correlation (such as the Gunnarsson-Lundqvist ||23| or Ceperly-Alder ||2^ prescrip- 
tions) to the effective HS potential results in an undesirable lowering of the peak-profile. If 
the correlation is added at right the level of exchange- correlation hole {p^ replaced by p^c in 
Eq. [§] ) and then the HS computations be performed self-consistently, an improvement over 
the "exchange only" J{q) is expected. On the other hand, though the work of Holas and 
March |T^ as pointed above prescribes an inclusion of the kinetic piece of correlation in the 
KS context it is not known how one actually implements their scheme in practice. These 
studies of course, go beyond the scope of the present "exchange-only" theme. 

In this paper, we have carried out a comparative study of the momentum space properties 
of atoms viz. Compton profiles and various expectation values calculated within the "work 
formalism" of Harbola and Sahni, the Hartree-Fock theory as well as the Self-Interaction- 
corrected Local Density Approximation theory. The Compton profiles for various excited 
states of the Helium atoms are also presented within the work formalism. The present work 
demonstrates that the Harbola-Sahni work formalism which in position space closely follows 
the HF theory also seen to do so in the momentum space. 
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TABLES 
TABLE L Spherically averaged Compton profile, Jsphil) for Helium within the three "exchange 

only " theories compared with their experimental counterpart. (Hartree a.u. used throughout) 



q 



HF 



HS 



SIC 



Expt."" 



0.0 
0.2 
0.6 
LO 
L5 
2.0 
2.5 
3.0 
"" Ref. M 



L070 


1.070 


1.070 


L017 


1.017 


1.017 


0.700 


0.700 


0.700 


0.382 


0.382 


0.382 


0.160 


0.160 


0.160 


0.068 


0.068 


0.068 


0.031 


0.031 


0.031 


0.015 


0.015 


0.015 



1.071±1.5% 
1.019 
0.705 
0.388 

0.069 

0.030±15% 

0.013 
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TABLE II. Jsphiq) for Neon 



q 


HF 


HS 


SIC 


CI" 


Expt.* 


0.0 


2.727 


2.719 


2.751 


2.739 


2.762 


0.2 


2.696 


2.687 


2.717 


2.707 


2.738 


0.4 


2.593 


2.585 


2.608 


2.602 


2.630 


0.6 


2.413 


2.406 


2.418 


2.4159 


2.427 


0.8 


2.168 


2.162 


2.163 


2.1645 


2.162 


1.0 


1.889 


1.885 


1.875 


1.880 


1.859 


1.5 


1.228 


1.228 


1.211 


- 


- 


2.0 


0.771 


0.774 


0.764 


0.768 


0.765 


2.5 


0.501 


0.506 


0.501 


- 


0.501 


3.0 


0.346 


0.350 


0.349 


0.348 


0.359 


3.5 


0.253 


0.256 


0.256 


- 


0.277 


4.0 


0.194 


0.196 


0.197 


0.196 


0.210 


5.0 


0.125 


0.125 


0.126 


0.126 


0.126 



" Ref . 0] 
'' Ref. [1191 
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TABLE III. Jsphiq) for Argon 



HF 



HS 



SIC 



Expt." 



0.0 


5.064 


0.2 


4.963 


0.4 


4.619 


0.6 


4.035 


0.8 


3.333 


1.0 


2.661 


1.5 


1.546 


2.0 


1.084 


2.5 


0.874 


3.0 


0.736 


3.5 


0.622 


4.0 


0.520 


4.5 


0.433 


5.0 


0.359 


10.0 


0.075 


15.0 


0.025 


'^Ref. 


15 



5.040 
4.941 
4.605 
4.029 
3.331 
2.664 
1.557 
1.090 
0.876 
0.736 
0.621 
0.519 
0.432 
0.359 
0.076 
0.025 



5.093 
4.991 
4.638 
4.033 
3.312 
2.636 
1.540 
1.086 
0.875 
0.737 
0.620 
0.520 
0.432 
0.359 
0.076 
0.025 



5.058 
4.917 
4.526 
3.960 
3.319 
2.697±1% 

1.129 
0.904 
0.744 
0.634 
0.534±2.5% 

0.366 

0.078±10% 
0.025 
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TABLE IV. Jsphiq) for Krypton 



q 


HF 


HS 


SIC 


Expt.° 


0.0 


7.237 


7.195 


7.262 


7.188 


0.2 


7.095 


7.060 


7.122 


6.988 


0.4 


6.605 


6.586 


6.625 


6.453 


0.6 


5.785 


5.781 


5.783 


5.702 


0.8 


4.855 


4.863 


4.841 


4.883 


1.0 


4.044 


4.059 


4.032 


4.131±1.7% 


2.0 


2.442 


2.447 


2.448 


2.557 


3.0 


1.858 


1.857 


1.854 


- 


4.0 


1.327 


1.324 


1.319 


1.350 


5.0 


0.935 


0.935 


0.931 


0.933±3.5% 


10.0 


0.260 


0.260 


0.261 


0.254 


15.0 


0.105 


0.105 


0.105 


0.099 


'^ Ref. 


15 









TAB 


LEV. 


< p" > moments for Helium within different "exchange-only" tlieories. 


moments 




HF HS SIC 



<p"^> 



<p ^> 



<p > 



<p^> 



<p^> 



<p'^> 



4.0893E+00 
2.1406E+00 
2.7990E+00 
5.7234E+00 
1.7991E+01 
1.0567E+02 



4.0760E+00 
2.1409E+00 
2.7990E+00 
5.7234E+00 
1.7990E+01 
1.0549E+02 



4.0902E+00 
2.1410E+00 
2.7987E+00 
5.7138E+00 
1.7628E+01 
8.7395E+02 
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TABLE VI. < p" > moments for Neon. 



moments 



HF 



HS 



SIC 



<p "■> 



<p ^> 



<p > 



<p^> 



<p^> 



<p'^> 



5.4795E+00 
5.4558E+00 
3.5196E+01 
2.5709E+02 
3.5843E+03 
9.8510E+04 



5.4526E+00 
5.4387E+00 
3.5269E+01 
2.5708E+02 
3.5720E+03 
9.9418E+04 



5.6349E+00 
5.5025E+00 
3.5246E+01 
2.5771E+02 
3.5836E+03 
9.9898E+04 



moments 



TABLE VII. The <p"> momets for Argon. 



HF 



HS 



SIC 



<p "■> 



<p ^> 



<p > 



<p'^> 



<p^> 



<p'^> 



1.3107E+01 
1.0128E+01 
8.8699E+01 
1.0536E+03 
2.4301E+04 
1.1393E+06 



1.2943E+01 
1.0076E+01 
8.8796E+01 
1.0536E+03 
2.4307E+04 
1.1723E+06 



1.3253E+01 
1.0187E+01 
8.8793E+01 
1.0538E+03 
2.3997E+04 
5.4391E+06 
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TABLE VIII. <p"> moments for Krypton. 



moments 



HF 



HS 



SIC 



<p "■> 



<p ^> 



<p> 



<p^> 



<p^> 



<p'^> 



1.7478E+01 
1.4474E+01 
2.8141E+02 
5.5041E+03 
2.2628E+05 
1.9852E+08 



1.7084E+01 
1.4390E+01 
2.8161E+02 
5.5013E+03 
2.2424E+05 
5.0212E+07 



1.7517E+01 
1.4524E+01 
2.8155E+02 
5.5072E+03 
2.2453E+05 
5.0297E+07 



TABLE IX. Jsphio) for different states of Helium atom within the Harbola-Sahni approach. 



q 


1^2 


ls2s 


ls2p 


ls3p 


ls4p 


0.0 


1.070 


2.516 


1.583 


2.966 


4.433 


0.2 


1.017 


1.532 


1.467 


1.266 


0.955 


0.4 


0.879 


0.592 


0.949 


0.517 


0.423 


0.6 


0.700 


0.362 


0.537 


0.412 


0.357 


0.8 


0.527 


0.294 


0.340 


0.303 


0.287 


1.0 


0.382 


0.237 


0.239 


0.228 


0.223 


1.5 


0.160 


0.119 


0.111 


0.112 


0.111 


2.0 


0.068 


0.056 


0.052 


0.053 


0.053 


2.5 


0.031 


0.029 


0.025 


0.025 


0.025 


3.0 


0.015 


0.013 


0.012 


0.012 


0.012 
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